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Introduction 


The  focus  of  this  research  is  to  determine  the  impact  of  our  knowledge 
of  ocean  variability,  as  characterized  by  the  sound  speed  profile,  on  ocean 
surveillance  detection  performance  for  situations  of  interest  to  the  U.  S.  Navy. 
Signal  detection  theory  provides  a  framework  for  determining  optimal  phase 
coherent  detection  and  source  localization  algorithms  and  determining  their 
performance.  This  framework  allows  us  to  incorporate  directly  the  physics  of 
acoustic  propagation  through  a  space-time  variable  medium  into  the  struc¬ 
ture  of  optimal  algorithms.  Thus  the  signal  processing  algorithms  designed 
from  this  viewpoint  have  a  built-in  robustness  to  the  environmental  vari¬ 
ability.  Using  these  optimal  algorithms  we  can  determine  the  maximum 
attainable  detection  performance. 

2  Approach 

Our  approach  is  a  combination  of  theory,  some  very  computationally  inten¬ 
sive  numerical  work,  and  the  use  of  real  sound  speed  profiles,  where  appropri¬ 
ate,  obtained  from  the  ocean  acoustics  community.  The  approach  is  to  use 
the  framework  of  signal  detection  theory  to  obtain  optimal  algorithms  for 
inhomogeneous  and  stochastic  ocean  models  characterized  by  an  ensemble  of 
sound  speed  profiles.  Ideally,  we  would  like  to  compute  the  detection  and  lo¬ 
calization  performance  for  a  volume  of  received  points  as  a  function  of  a  large 
volume  of  possible  source  locations,  for  varying  degrees  of  knowledge  of  an 
uncertain,  inhomogeneous  and  stochastic  ocean  medium.  Available  software 
codes  and  computational  speeds  limit  the  complexity  of  situations  that  can 
be  considered.  The  approach  requires:  (a)  setting  up  a  meaningful  sound 
speed  profile  description  of  the  ocean;  (b)  incorporating  the  variabilities  and 
stochastic  nature  of  the  ocean  in  the  description  of  the  ensemble  of  sound 
speed  profiles,  as  well  cis  uncertainties  of  source  signal  characteristics,  source 
position,  and  noise;  (c)  solving  for  the  received  fields  as  a  function  of  all  pos¬ 
sible  source  locations;  (d)  using  the  received  fields  to  generate  the  optimal 
detection  and  localization  algorithms  for  a  given  ocean  acoustic  scenario  (i.e. 
the  likelihood  ratio  for  the  detection  problem  and  the  a  posteriori  probability 
for  the  source  localization  and  the  tomography  problem);  and  (e)  obtaining 
the  detection  and  localization  performance  of  these  optimal  algorithms  to  get 
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the  maximum  obtainable  performance. 

Our  approach  has  been  to  exploit  the  ability  to  remotely  sense  the  sound 
velocity  profile  of  the  ocean,  as  demonstrated  by  the  ocean  acoustic  tomogra¬ 
phy  group  of  Scripps/Woods  Hole/Univ.  of  Michigan/APL.  A  very  important 
issue  is  to  determine  what  the  attainable  surveillance  detection  performance 
is  given  that  one  has  a  more  detailed  description  of  how  sound  is  propagated 
in  an  inhomogeneous  ocean  environment,  in  the  deep  ocean,  at  long  ranges. 
Signal  detection  theory  provides  a  framework  for  determining  optimal  detec¬ 
tion  algorithms  and  their  performance.  In  other  words,  we  are  able  to  utilize 
the  knowledge  that  the  ocean  tomography  group  is  able  to  obtain  in  regard 
to  the  sound  speed  profile  and  use  that  physics  to  determine  the  changing 
wavefront  of  the  propagating  signal  in  an  inhomogeneous  ocean  environment. 
This  information  is  used  to  drive  the  structure  of  optimal  algorithms  from 
which  we  can  obtain  the  maximum  attainable  detection  performance. 


3  Summary  of  Research  Accomplished 

Results  have  been  obtained  in: 

•  Development  of  new  a  posteriori  probability  source  localization  algo¬ 
rithms  by  Richardson  and  Nolte  [10,  11].  These  algorithms,  unlike 
matched  field  processing  algorithms,  do  not  require  complete  knowledge 
of  the  acoustic  environment,  but  can  determine  source  position  even 
with  uncertain  or  imprecise  information  about  the  environment.  This 
algorithm  is  termed  the  optimum  uncertain  field  processing  (OUFP) 
algorithm. 

Paramer  estimation  theory  is  utilized  to  derive  the  new  algorithm.  THis 
provides  a  systematic,  optimal  approach  to  the  problem,  and  allows  en¬ 
vironmental  uncertainty  to  be  easilyincorporated  into  the  algorithm.  In 
addition  to  estimating  source  position,  estimates  of  parameters  of  the 
acoustic  environment  can  also  be  calculated.  This  makes  simultane¬ 
ous  source  localization  and  acoustic  tomographic  estimation  of  ocean 
parameters  possible. 

A  series  of  results,  some  of  which  are  presented  here,  illustrate  the 
robust  performance  of  the  uncerain  field  processor,  relative  to  the  per- 
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formance  of  matched  field  processing  methods.  Estimation  of  ocean 
acoustic  parameters  is  also  illustrated. 

•  Development  of  optimum  ocean-model-based  detection  algorithms  by 
LazofF  [8].  Effects  of  signal  bandwidth  on  detection  performance  of 
optimal  ocean-model-based  algorithms,  for  signals  received  through  an 
inhomogeneous  acoustic  ocean  medium. 

•  Sensitivity  of  wide  band  matched-field  localization  to  the  sound  speed 
profile  information  by  Clarke  and  Rausch  [2,  12]. 
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4  Results-Optimal  Acoustic  Source  Localiza¬ 
tion  Algorithm 

The  problem  of  localization  in  cases  in  which  there  is  uncertainty  about  the 
parameters  describing  the  acoustic  environment  is  considered  by  Richard¬ 
son  [10,  11].  This  is  a  generalization  of  the  problem  in  which  it  is  desired 
to  locate  a  source  in  a  known  or  certain  environment  (the  matched  field 
problem),  which  has  received  considerable  attention  in  the  literature  since 
Bucker’s  1976  paper  [Ij. 

The  sensitivity  of  matched  field  algorithms  to  environmental  uncertainty 
has  been  illustrated  in  various  studies  [1,  9,  3,  13,  10].  These  studies  in¬ 
dicate  that  matched  field  localization  algorithms  can  be  quite  sensitive  to 
mismatch  between  the  actual  environment  and  that  for  which  the  algorithms 
are  designed. 

By  considering  the  more  general  problem,  we  can  design  new  algorithms 
which  we  hope  are  less  sensitive  to  mismatch.  Since  the  known  environment 
is  a  limiting  case  of  the  uncertain  environment,  valuable  insight  is  gained 
into  the  matched  field  problem. 

The  source  localization  problem  is  regarded  as  a  problem  in  parameter  es¬ 
timation  theory.  The  source  position  coordinates  (range,  depth  and  bearing) 
are  treated  as  unknown  parameters  of  the  observed  acoustic  field.  Parame¬ 
ter  estimation  theory  provides  the  necessary  framework  for  the  design  of  an 
optimum  algorithm  for  estimating  source  position.  This  algorithm  is  termed 
the  optimum  uncertain  field  processor  (OUFP).  The  main  advantages  of  this 
approach  are;  (1)  it  provides  a  framework  for  the  development  of  optimum 
localization  algorithms,  (2)  it  introduces  well  understood  performance  mea¬ 
sures  for  comparison  of  different  algorithms,  and  (3)  it  allows  probabilistic 
descriptions  of  environmental  uncertainty  to  be  easily  incorporated  into  the 
problem. 

4.1  Notation 

The  observation  (or  measurement)  consists  of  K  time  samples  from  each  of 
N  elements  of  a  receive  array.  The  vector  Vk  represents  a  snapshot  of  the 
element  outputs  at  time  sample  k.  The  complete  observation  is  arranged  as 
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an  NK  X  1  matrix  (column  vector), 

r  =  [r^rj; . . .  !r]:r  (1) 

The  observation  consists  of  signal  and  additive  noise  components, 

r  =  s(S,  'i')  +  w.  (2) 

We  assume  that  the  signal  component  depends  upon  several  parameters, 
some  of  which  are  unknown.  The  set  of  source  position  parameters  (e.g.. 
range,  depth,  and  bearing)  is  denoted  S.  Any  parameters  describing  the 
transmitted  pressure  waveform  (e.g.,  amplitude,  phase,  center  frequency, 
etc.)  are  denoted  Parameters  which  are  required  to  describe  the  acous¬ 
tic  medium  (e.g.,  any  sound  velocity  profile  parameters,  descriptions  of  the 
boundary  conditions,  etc.)  are  collectively  denoted  'i'. 

4.2  The  General  OUFP  Equation 

The  range,  depth  and  bearing  of  the  source  are  treated  as  random  variables. 
The  optimum  uncertain  field  processor  is  defined  as  the  implementation  of 
the  source  position  a  posteriori  probability  density  function.  Source  and  envi¬ 
ronment  uncertainty,  as  expressed  by  probability  density  functions,  can  easily 
be  included  into  the  problem.  This  definition  is  intuitively  meaningful  and 
allows  many  of  the  classical  estimators  (for  example,  minimum  mean-square 
error,  maximum  a  posteriori,  maximum  likelihood,  etc.)  to  be  calculated. 

The  source  position  a  posteriori  probability  density  function  (pdf),  ps|r(S|r), 
is  related  to  the  conditional  observation  pdf,  pr|s(r|S)  by 

rei,(S|r)  =  ^p,|s(r|S).  (3) 

Prir) 

Only  the  case  in  which  the  observation  is  a  deterministic  function  of  the 
transmitted  waveform  and  S  and  'i'  is  considered  here  (the  deterministic 
ocean  model).  The  transmitted  signal  is  assumed  to  be  completely  deter¬ 
mined  if  #  is  known  (the  deterministic  signai  ca.se).  The  pdf  describing  the 
observation  is  then  related  to  the  noise  pdf  by 

Pr|S,^,'i'(rjS,  #,  4^)  =  pw{r  -  s(S,  $,  '^)).  (4) 
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(This  e  , nation  is  not  applicable  if  either  the  transmitted  signal  or  the  acoustic 
transfer  function  are  random). 

From  the  definition  of  conditional  probability  density  functions  we  can 
write 

Pr,$,'®is(r,  'i'IS)  =  Pris.*,'*'(r|s,  ^|S),  (5) 

or 

Pr,'*,*|s(r,  ^|S)  =  Pw(r  -  s(S,  '®))p$,^-|s($,  ’®'|S).  (6) 

We  can  obtain  pr|s(r|S)  by  integrating  this  density  function  over  $  and 

Pr|s(r|S)  =  Pw(r  -  s(S,  $,  ^))p$,^is($,  ^|S)  d^.  (7) 

Finally,  the  desired  a  posteriori  source  position  pdf  can  be  obtained  by 
substitution  of  this  expression  into  Eq.  3, 

ra|r(S|r)  =  ^  y^/^p„(r-s{S,#,®))p*,»|s(«,®|S)d$rf».  (8) 

This  is  the  defining  equation  of  the  optimum  uncertain  field  processor. 
It  can  be  simplified  in  many  applications  of  interest,  as  shown  in  the  next 
subsection. 

4.3  The  Narrowband  Signal  in  Gaussian  Noise  Prob¬ 
lem 

Matched  field  processing  methods  have  concentrated  on  localizing  a  nar¬ 
rowband  source  in  the  presence  of  a  Gaussian  noise  background.  It  is  this 
problem  to  which  we  now  turn  our  attention. 

If  the  elements  of  the  noise  component  of  the  observation  are  Gaussian 
distributed  and  jointly  stationary,  the  noise  probability  density  function  is 
given  by 

fl  is  the  space-time  noise  covariance  matrix  and  |f2|  is  its  determinant.  The 
conditional  pdf  of  the  observation  can  therefore  be  written  as 

exp  (— 

Pr|S.*.«'(r|S,#,^)  = 
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where 


A(S,$,'®',r)  =  is^O-‘s.  (11) 

£t 

(For  convenience,  the  explicit  notational  dependence  of  s  upon  S,  and 
has  been  dropped.) 

The  source  position  a  posteriori  pdf  is  equal  to 
PS|r(Slr)  =  C(r)ps(S)  exp  (A(S,  $, r))  |S)  d$ 

(12) 

where  C(r)  is  some  function  of  r  (the  functional  relationshio  is  of  no  signif¬ 
icance)  that  reduces  to  a  normalization  constant  for  any  particular  observa¬ 
tion. 

The  problem  is  simplified  if  the  noise  samples  are  uncorrelated  in  time. 
The  matrix  fZ  can  then  be  written  as  the  Kronecker  product  of  the  spatial 
covariance  matrix,  Q,  and  a  K  x  K  identity  matrix,  IkKi  that  is, 

fZ  =  Q  (13) 

(The  Kronecker  product  is  the  NO  X  MP  matrix  which  results  from  the 
element  by  element  multiplication  of  an  iV  x  M  by  an  O  x  P  matrix.)  The 
following  set  of  equations  are  easily  derived  from  the  properties  the  Kronecker 
product. 


fZ-‘  = 

Q-’  ®  Ikk 

K 

(14) 

r’’n-‘s  = 

k=l 

K 

(15) 

s^fZ~'s  = 

fc=l 

(16) 

Eq.  11  can  then  be  written  as 

A(S,$,^,r)  =  f: 

Jt=i 

’sfc  -  ^s[Q  ^Sk 

(17) 

An  eigenvector  decomposition  of  the  spatial  covariance  matrix  can  be 
represented  by 

N 

Q  =  S  ^nqnq;, 

n=l 

(18) 
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where  the  q„  are  the  eigenvectors  of  Q  and  the  A„  are  the  corresponding 
eigenvalues.  This  can  be  used  to  rewrite  Eq.  17  as 

A(S,$,’®f,r)  =  ^  A-*q^  H  (sfcr[  -  q„.  (19) 

n=l  U=1  ^  ^  . 

A  transmitted  narrowband  pressure  waveform  has  the  form 

x{i)  =  Acos(a;oi  +  &)■  (20) 

The  transmitted  amplitude  and  phase  are  assumed  to  be  unknown  parame¬ 
ters  of  the  transmitted  waveform;  i.e.,  they  are  members  of  the  set  $. 

The  acoustic  transfer  function  between  the  source  and  array  element  n 
can  be  written  as 

//„(S,^)  =  a„(S,'®)e‘®'‘<S-«-).  (21) 

The  transfer  function  depends  only  upon  the  source  position  and  the  param¬ 
eters  of  the  environment. 

The  signal  component  of  the  observation  is  therefore  equal  to 


Snk  =  AOn  COs(u;o4  -|-  +  &)■  (22) 

The  outer  product  terms  in  Eq.  19  can  now  be  evaluated.  The  element 
of  the  matrix  in  row  n  and  column  n'  is  equal  to 

■  K'  1  K 

y~!  =  A^Gnan'  COs(u;of;t  +  +  ^)  COs(u;oiit  +  On’  +  0).  (23) 

•^—1  ■  nn'  /c=l 

If  the  sampling  interval  is  much  less  than  the  signal  period,  the  summation 
may  be  approximated  as  an  integral.  This  results  in 


- - J  cos(u;ot  O-n  0)  cos(u;ot  -j-  On'  4"  0)  dt,  (24) 


where  T  is  the  observation  time.  A  trigonometric  identity  reduces  this  to 


A^KanOn’ 

2r  ^ 


—  On')  4"  cos(2u7ot  -f-  20  4"  On  4"  On')  dt. 


(25) 
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Assuming  that  the  observation  time  is  over  several  periods  of  the  observation, 
the  integral  over  the  second  term  is  negligible  when  compared  to  the  first  and 
we  have 

jr  (IjiOtiI  COs(^0ji  ^n')‘  ("fi) 


*=1 


Rewriting  this  in  terms  of  the  acoustic  transfer  function  yields 

A^K, 


r  K 

k=l 


The  entire  matrix  can  therefore  be  written  as 

fc=l  ^ 


(27) 


(28) 


where  H  is  an  element  column  vector  whose  elements  are  given  by  Eq.  21. 

Turning  to  the  other  outer  product  term,  an  element  of  the  matrix, 
HfcLi  is  equal  to 


K 

.*=1 


K 


=  ^  H  f’n>{tk)an  COs{u!otk  +  On  +  9). 

k=l 


Again  approximating  the  summation  by  an  integral  results  in 


.k=i 


^  fT' 
T 


f  r„.(t)a„  cosfwo^  +  On  +  0)  dt, 

Jo 


(29) 


(30) 


which  can  be  written  in  complex  notation  as 


r  K 


U=i 


=  (31) 


The  notation  can  be  simplified  by  rewriting  this  as 


r  K 

k=l 


A(TaK 


where  Pn'  is  defined  as 


Pn'  = 


2^  F 

(TaT 


I  rn'{t)t 

Jo 


-.wot 


(32) 


(33) 
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and  (Ta  is  a  normalization  factor  that  will  be  defined  shortly.  The 
matrix  is  therefore  equal  to 

complete 

(34) 

The  vector  P  is  an  TV  element  column  vector  whose  elements  are  defined  in 
Eq.  33. 

Substituting  Eqs.  28  and  34  into  Eq.  19  and  recalling  Eq.  18  results  in 

A(S, «,  ®,  r)  =  -^E  +  ^n{Re“} , 

(35) 

where 

E  = 

(36) 

R  = 

(37) 

R  indicates  to  what  degree  the  measured  field,  P,  is  similar  to  the  matching 
field,  H.  £  is  a  normalization  term  which  is  proportional  to  the  self-energy 
in  the  matching  field.  Both  R  and  E  are  functions  of  S,  $  and  'fp  through 
their  dependence  on  H.  If  we  define  a  detection  factor,  F  as 

b 

II 

(38) 

then  E  and  R  can  be  rewritten  as 

E  = 

R  =  cr^FH-Q-^P. 

(39) 

(40) 

If  R  is  expressed  in  polar  notation  as  /?  =  Eq.  35  can  be  written 

in  terms  of  real  variables  as 

A(S,  #,  r)  =  -^E  -b  — |F|  cos{e  +  Or) 

Z(T%  <Ja 

(41) 

It  is  assumed  that  #  and  ^  are  statistically  independent;  the  values  of 
amplitude  and  ph<ise  of  the  acoustic  source  have  no  dependence  on  parame¬ 
ters  describing  the  environment.  Furthermore,  the  amplitude  is  statistically 
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independent  of  the  phcise,  and  both  amplitude  and  phase  are  statistically 
independent  of  source  position.  These  assumptions  allow  the  probability 


density  function  in  Eq.  12  to  be  written  as 

P^,'*'|s(^,^|S)  =  p>»(A)p0(^)p$|s(’i^|S).  (42) 

Thus  the  OUFP  equation  becomes 

ps|r(S|r)  =  C(r)ps(S)  n(S,  r)p^|s(’^|S)  (43) 

where 

n(S,^,r)=  /  f  expmS,^,d,A,T))pA{A)pe{9)d9dA  (44) 

JA  Je 

Maximum  uncertainty  is  assumed  about  the  transmitted  phase,  i.e.,  it  is 
described  by  the  following  probability  density  function, 

PtW  =  ii  0<e<2T.  (45) 

From  reference  [7,  Eq.  3.915.4]  we  have 

exp(/?cos(a:  +  0))  dx  = /o(/9),  (46) 

Itt  Jo 

where  Iq{I3)  is  a  modified  Bessel  function  of  order  zero.  Using  this  equation 
and  Eqs.  41  and  45  in  Eq.  44  gives 

n(S,»,r)=  PA{A)dA.  (47) 

If  the  amplitude  is  assumed  to  be  Rayleigh  distributed  with  pdf 

Pa{A)  =  ^exp(^-^^  A>0,  (48) 

a  closed  form  expression  for  Eq.  47  is  given  in  reference  [4,  page  311]: 

n(S,®,r)  =  (£:+l)-'expf^j.  (49) 
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Collecting  everything  together,  the  OUFP  equation  for  the  narrowband 
signal  in  Gaussian  noise  problem  is 

ps|r(S|r)  =  C(r)ps(S)  j^{E  +  1)-'  exp  pt|s(*|S)d«,  (50) 

where 

E  =  (tIFWQ-^H  (51) 

R  =  <tIFWQ-^P.  (52) 

4.4  Matched  Field  Processing 

If  the  acoustic  environmental  parameters  are  known,  the  optimum  uncertain 
field  processor  reduces  to 

rei,(S|r)  =  C(r)ps(S)(E  +  1)-'  exp  ,  (53) 

where  E  and  R  are  now  calculated  only  for  the  known  environment.  This 
is  the  optimum  algorithm  for  locating  a  narrowband  source  in  a  known  en¬ 
vironment  and  the  implementation  of  this  equation  is  termed  the  optimum 
matched  field  processor  (OMFP). 

It  is  instructive  to  consider  maximum  a  posteriori  (MAP)  estimation  of 
the  source  position.  MAP  estimates  are  obtained  by  determining  the  peak  of 
pS|r(S|r).  Equivalently,  we  can  determine  the  peak  of  any  monotonic  function 
of  ps|r(S|r).  A  particularly  illuminating  monotonic  function  is 

/(S)  =  i  In  (ps|,(S|r)/C(r))  .  (54) 

Substituting  from  Eq.  53  gives 

/(S)  =  jJ|!^-i|n(E+l)  +  );ln(ps(S)).  (55) 

Using  the  definitions  of  R  and  E,  and  assuming  that  E  >•  1,  yields 

/(S)  =  +  j  ln(ps(S)).  (56) 
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For  high  values  of  the  detection  factor  this  can  be  approximated  as 


.2  1h-q-^p|^ 

“  H’Q-IH 

and  if  the  noise  is  spatially  uncorrelated  this  reduces  further  to 


(57) 


IH*PP 

This  equation  was  first  suggested  for  use  in  matched  field  localization  in 
reference  [6].  It  is  a  generalization  of  the  conventional  (Bartlett)  beamformer. 
(In  conventional  beamforming  H  corresponds  to  the  steering  vector  for  a 
plane- wave  arrival.)  The  surface  obtained  by  plotting  /(S)  as  a  function 
of  S  for  a  particular  (usually  noise  free)  observation,  P,  is  known  as  an 
ambiguity  surface. 


4.5  Acoustic  Tomography 

The  source  position  a  posteriori  probability  density  function  can  also  be 
written  as 

PSir(S|r)  =  PS,'® (S,  ^\t)  d^.  (59) 

That  is,  the  OUFP  equation  is  a  marginal  density  function  computed  from 
the  joint  probability  density  function  of  S  (the  set  of  source  position  coordi¬ 
nates)  and  ^  (the  set  of  environment  parameters). 

From  equation  50  the  joint  pdf  can  be  identified  as 

PS.«(S,  *|r)  =  C{r)ps(S)pi{^)(B  +  1)-'  exp  ■  (60) 

(S  and  $  are  cissumed  to  be  statistically  independent.)  ps(S)  and  p®(^)  are 
the  a  priori  probability  density  functions.  Recall  that  E  and  R  are  dependent 
on  S  and  through  their  dependence  on  the  acoustic  transfer  function,  H. 

It  has  been  shown  that,  if  ^  is  known,  this  equation  reduces  to  the 
optimum  matched  field  equation.  If  there  is  uncertainty  in  then  computing 
the  marginal  density  over  ^  results  in  the  optimum  uncertain  field  processor. 
Other  situations  are  worthy  of  mention.  If  the  source  position  is  known, 
or  the  marginal  over  the  source  position  is  calculated,  then  the  join*^  pdf 
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reduces  to  an  a  posteriori  pdf  describing  the  environment  parameters, 

The  pdf  can  be  used  to  estimate  the  values  of  the  environment  parameters. 
From  these  estimates  a  “picture”  of  the  acoustic  sound  speed  field  can  be 
obtained.  This  represents  one  approach  to  the  acoustic  ocean  tomography 
problem.  Another  concept  suggested  by  the  previous  equation  is  that  of 
simultaneous  localization  and  tomography.  Some  examples  of  this  are  given 
in  the  results. 
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5  Results-Sensitivity  to  Large  Scale  Range 
Independent  Perturbations 

In  this  section,  localization  processor  sensitivity  to  uncertainty  in  the  ma¬ 
jor  parameters  of  the  deep-ocean  sound  velocity  profile  is  investigated  by 
Richardson  [10,  11].  A  subsection  on  acoustic  tomographic  estimation  of 
sound  velocity  profile  parameters  is  also  included. 

A  range-independent  acoustic  environment  is  simulated  and  the  normal 
mode  acoustic  propagation  model  is  used. 

5.1  The  Sound  Velocity  Profile  Model 

The  sound  velocity  profile  is  modeled  as 

c(z)  =  cm(2)  +  wiei(z) -f  1^262(2),  (61) 

where  cm{z)  is  a  Munk  profile  with  parameters  cq  (axis  sound  speed)  equal  to 
1500  m/s  and  zq  (axis  depth)  equal  to  1000  m.  ei{z)  and  e2{z)  are  orthogonal 
perturbation  functions.  We  will  investigate  the  performance  sensitivity  of  two 
localization  processors  with  respect  to  uncertainty  in  the  weight  parameters 
lUi  and  W2. 

The  functions  e\{z)  and  62(2)  were  calculated  from  a  Giam-Schmidt  or- 
thogonalization  [5]  on  the  function  set  {/^(z), 72(2)},  where 

h{z)  =  ci{z)-cm{z)  (62) 

f2{z)  =  C2(z)  -  cm{z).  (63) 

ci(z)  and  C2{z)  are  Munk  profiles  with  parameters,  cq  =  1510  m/s  (axis 
speed),  Zq  =  1000  m  (axis  depth)  and  Cq  =  1500  m/s,  zq  =  950  m  respectively. 
Cl (2)  differs  from  cm{z)  only  in  the  value  of  the  cq  parameter,  while  02(2) 
differs  from  c\f{z)  only  in  the  value  of  the  zq  parameter. 

The  profile  cm{z)  is  plotted  in  the  center  graph  in  Fig.  1.  The  other 
graphs  in  the  figure  are  plots  of  A{z)  =  wi€i{z)  +  W2€2{z)  for  different  values 
of  wi  and  u>2.  A  value  of  wi  other  than  zero  essentially  adds  a  constant 
value  to  cm{z)  for  all  depths.  Values  of  W2  different  from  zero  affect  only  the 
portion  of  the  SVP  above  the  channel  axis  and  can  be  used  to  raise  or  lower 
the  depth  of  the  sound  speed  axis. 
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5.2  Estimation  Sensitivity 

The  performance  of  the  optimum  uncertain  field  processor  (OUFP)  and  that 
of  the  optimum  matched  field  processor  (OMFP)  are  examined  in  the  pres¬ 
ence  of  sound  velocity  profile  uncertainty. 

Since  the  optimum  matched  field  processor  is  not  the  optimum  localiza¬ 
tion  processor  to  use  in  an  uncertain  environment,  the  terminology  is  a  little 
confusing.  The  word  “optimum",  in  the  label  “optimum  matched  field  pro¬ 
cessor”,  is  now  being  used  to  identify  a  particular  matched  field  processor 
and  not  as  an  indication  of  level  of  performance.  To  remedy  the  problem,  the 
notation  mismatched-OMFP  will  be  used  to  indicate  the  optimum  matched 
field  processor  when  operating  in  an  uncertain  environment. 

5.2.1  Range/Depth  Estimation 

In  the  simulated  results  presented  in  this  section  the  acoustic  source  is  at 
a  depth  of  1000  m  and  at  a  range  of  50000  m  from  the  receive  array.  A 
vertical  array  of  10  elements,  with  the  shallowest  element  at  200  m  and  the 
elements  spaced  200  m  apart,  is  used.  The  transmitted  frequency  is  15  Hz. 
The  detection  factor  is  92  dB. 

The  mismatched-OMFP  calculates  replica  fields  assuming  that  and 
W2  are  both  zero.  (It  is  matched  to  only  the  environment  characterized  by 
those  values.) 

The  OUFP  assumes  that  w\  is  uniformly  distributed  between  the  values 
-0.2  and  0.2  and  that  tnj  is  equal  to  zero. 

We  investigate  processor  performance  by  computing  and  plotting  the  a 
posteriori  source  position  probabilities.  Each  point  in  the  a  posteriori  prob¬ 
ability  graphs  represents  the  probability  that  the  acoustic  source  is  located 
within  a  5  meter  depth  bin  and  a  1000  meter  range  bin. 

The  output  of  the  mismatched-OMFP  for  an  actual  environment  that  is 
characterized  by  values  of  wi  and  W2  both  equal  to  zero  is  shown  is  shown  in 
Fig.  2.  This  is  the  environment  for  which  the  processor  is  designed.  It  does 
a  good  job  of  locating  the  source  as  evidenced  by  the  small  cluster  of  peaks 
near  the  true  source  position. 

The  output  of  the  OMFP  for  the  same  environment  is  shown  in  Figure  3. 
The  observation  (signal  and  noise)  is  the  same  as  that  for  the  mismatched- 
OMFP  of  the  previous  figure.  The  two  processors  do,  of  course,  process 
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the  observation  differently.  The  OUFP  also  does  a  good  job  of  locating  the 
source;  however,  the  main  peak  is  slightly  lower  than  that  in  the  previous 
figure  and  some  small  spurious  peaks  at  a  distance  away  from  the  source. 
This  can  be  explained  by  the  fact  that  the  OUFP  is  looking  over  a  larger 
set  of  possible  environments  than  the  mismatched-OMFP  and  there  is  a 
corresponding  penalty  in  performance.  The  OUFP  may  not  be  the  best 
processor  in  any  particular  environment,  but  rather,  it  is  the  best  processor 
over  all  possible  uncertain  environments. 

The  output  of  the  mismatched-OMFP  for  an  actual  environment  char¬ 
acterized  by  w\  equal  to  0.2  and  W2  equal  to  0.0  is  shown  in  Fig.  4.  (This 
corresponds  to  a  mismatch  between  the  assumed  and  actual  sound  velocity 
profiles  of  approximately  2  m/s.)  The  processor  is  the  same  as  that  used  in 
Fig.  2,  but  the  observation  is  now  different.  The  processor  is  mismatched  to 
this  environment  and  the  localization  performance  is  poor.  There  is  no  peak 
at  the  actual  source  position,  although  some  small  peaks  appear  nearby,  and 
the  largest  peak  is  at  a  range  of  about  30000  m  and  a  depth  of  800  m. 

The  output  of  the  OUFP  for  the  same  observation  is  shown  in  Fig.  5. 
The  OUFP  does  a  much  better  job  of  locating  the  source  as  indicated  by  the 
large  peak  at  the  true  source  position.  (The  fact  that  performance  is  better 
here  than  in  Fig.  3  violates  no  principles.  The  processors  are  identical,  but 
the  observations  are  different.)  This  clearly  illustrates  the  advantage  of  using 
an  uncertain  field  processor  in  an  uncertain  environment. 

5.2.2  Range  Estimation 

We  now  concentrate  on  the  estimation  of  only  one  source  position  coordi¬ 
nate  at  a  time.  This  allows  us  to  examine  the  sensitivity  of  the  localization 
processors  to  uncertainty  in  the  environment  in  greater  detail. 

The  range  estimation  problem  is  considered  first.  The  source  depth  is 
known  and  equal  to  1000  m.  The  localization  scenario  is  the  same  «is  in  the 
previous  section,  except  that  the  detection  factor  has  been  reduced  to  86  dB. 
(The  higher  detection  factor  is  not  necessary  for  estimation  of  a  single  source 
coordinate.)  These  simulations  are  noise-free.  This  prevents  misinterpreta¬ 
tion  of  spurious  noise  peaks.  This  could  also  be  accomplished  by  looking  at 
many  noisy  simulation  runs,  but  this  is  prohibited  by  computer  processing 
time  requirements. 

The  localization  processor  output  is  now  an  a  posteriori  probability  line 
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instead  of  a  surface  as  in  the  previous  section.  Multiple  lines,  corresponding 
to  different  actual  environments,  are  displayed  simultaneously  in  a  three- 
dimensional  graph.  This  allows  us  to  visualize  the  change  in  processor  per¬ 
formance  with  a  change  in  environment.  In  an  actual  implementation  of 
the  localization  processor  only  one  line  (the  line  corresponding  to  the  actual 
environment)  will  be  observed. 

The  range  a  posteriori  pdf  output  of  the  mismatched-OMFP  for  different 
environment  values  of  lOj  is  shown  in  Fig.  6.  The  value  of  w\  is  zero.  This 
processor  is  matched  to  an  environment  characterized  by  values  of  w\  and 
W2  equal  to  zero.  The  figure  indicates  that  the  performance  of  the  processor 
is  not  sensitive  to  variation  in  between  -1  and  1.  (This  corresponds  to  a 
variation  in  the  depth  of  the  sound  axis  of  about  100  meters.) 

The  corresponding  result  for  the  OUFP  is  shown  in  Fig.  7.  The  OUFP 
assumes  that  the  value  of  W2  is  uniformly  distributed  between  —1  and  1. 
This  processor  is  also  not  sensitive  to  the  actual  value  of  'W2. 

The  mismatched-OMFP  output  for  different  values  of  the  environment 
parameter  loi  is  shown  in  Fig.  8.  The  value  of  W2  is  zero.  The  processor  is 
matched  to  an  environment  characterized  by  values  of  wi  and  W2  equal  to 
zero  and  there  is  a  large  peak  at  the  true  source  range  for  the  matched  case 
{wi  equal  to  zero).  In  the  mismatched  cases  however,  several  spurious  peaks 
start  to  appear  at  large  distances  from  the  true  source  position.  There  are 
cases  in  which  mismatch  is  less  than  .1  (corresponding  to  knowledge  of  the 
SVP  to  within  1  m/s)  and  the  largest  peak  occurs  at  a  distance  of  greater 
than  10000  meters  from  the  true  source  range. 

The  reduced  sensitivity  of  the  OUFP  is  illustrated  in  Figure  9.  This 
processor  assumes  that  the  value  of  lUx  is  uniformly  distributed  between  —1 
and  1.  A  distinctive  ridge,  extending  across  all  values  of  wi,  is  present  at 
the  true  source  range  of  50000  m.  Such  ridges  are  characteristic  of  robust 
localization  processors. 

5.2.3  Depth  Estimation 

We  now  consider  the  problem  of  estimating  the  acoustic  source  depth  when 
the  range  to  the  source  (50000  m)  is  known. 

The  output  of  the  mismatched-OMFP  as  a  function  of  the  parameter 
W2  is  shown  in  Fig.  10.  The  mismatched-OMFP  assumes  the  environment 
is  known  and  is  characterized  by  values  of  wi  and  W2  equal  to  zero.  The 
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Figure  7:  OUFP  a  posteriori  range  pdf  output  as  a  function  of  the  value  of 
102 .  wi  is  equal  to  zero. 
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result  indicates  that  depth  estimation  by  the  mismatched-OMFP,  like  range 
estimation,  is  not  sensitive  to  uncertainty  in  W2. 

The  corresponding  output  of  the  optimum  uncertain  field  processor  is 
shown  in  Fig.  11.  Depth  estimation  by  the  OUFP  is  also  not  sensitive  to 
uncertainty  in  w^.  The  height  of  the  ridge  is  just  slightly  less  than  that  in 
the  previous  figure. 

The  output  of  the  mismatched-OMFP  for  different  values  of  wi  is  shown 
in  Fig.  12.  There  is  a  cluster  of  peaks  at  the  true  source  depth  in  the  region 
where  loi  is  equal  to  zero.  For  values  of  wi  away  from  zero  spurious  peaks 
appear  at  depths  different  than  the  true  depth. 

The  OUFP  output  is  shown  in  Fig.  13.  A  ridge  at  the  true  source  depth, 
indicating  the  robustness  of  the  OUFP  when  operating  in  an  uncertain  en¬ 
vironment,  is  clearly  observed. 


5.3  Acoustic  Tomography 

As  discussed  in  Chapter  4,  a  processor  which  can  simultaneously  be  used 
to  locate  a  source  and  estimate  the  acoustic  environment  parameters  can 
be  constructed  based  on  Eq.  60.  Simulation  results  for  this  problem  are 
presented  in  this  section. 

The  source  position  and  receiver  geometry  are  the  same  as  in  the  previous 
section.  The  source  acoustic  frequency  is  15  Hz.  The  source  detection  factor 
is  86  dB. 

Simultaneous  estimation  of  one  source  position  coordinate  and  the  envi¬ 
ronment  parameter  is  considered  first.  The  other  position  coordinate  and 
W2  are  known. 

The  localization/tomography  processor  output  when  the  source  range  and 
the  value  of  w\  are  unknown  is  shown  in  Fig.  14.  (Recall  that  the  processor 
output  is  the  joint  range  and  wx  pdf.)  The  graph  indicates  that  there  is  a 
high  probability  that  the  source  is  located  at  a  range  of  50000  m  and  that 
the  value  wx  is  equal  to  zero.  Since  this  corresponds  to  the  true  values  of 
these  parameters,  it  can  be  concluded  that  the  processor  h«is  done  a  good 
job  in  this  caise. 

The  joint  depth  and  wx  pdf  is  shown  in  Fig.  15.  The  processor  also  does 
a  good  job  of  estimating  both  source  depth  and  lUi.  The  figure  indicates  a 
that  there  is  a  high  probability  of  locating  the  source  at  a  depth  of  1000  m 


29 


Loc.  Prob. 


ooo 

Depth 


of  tbe  value  of  uJ2.  en 


Depth 


Figure  U 
102.  ®' 


•  OUFP  a  posteriori  depth  pdf  outpu 


function  of  the  value  of 


qual  to 


zero. 


30 


Loc.  Prob. 


Deptln 


,  ,  omFP  a  ” 

Figure  12-.  So 

Xe  value  o!»..u.r»  equal  to  «re. 


g 

Q- 


1  ooo 

Oeptin 


OXJFP  a  posteriori  depth  pdf  output  as 


function  of  the  value  of 


Figure  13 
u,,.  W2  is  equal  to  zero. 


31 


Loc.  Prob. 


20 


A-'O  60 

Range 


depth  and  «  are  toow.. 


,5.  Localiaation/tomography  processor  joint  depth  and  .r  pdf  out 

The  range  and  uij  are  known. 


32 


and  that  the  value  of  wi  is  zero.  These  are  the  correct  values  of  these  two 
parameters. 

These  results  help  to  explain  why  the  mismatched-OMFP  is  sensitive  to 
uncertainty  in  wi  while  the  OUFP  is  not.  Considering  only  range  estimation 
for  a  moment  and  referring  to  Fig.  14,  the  output  of  the  mismatched-OMFP 
will  be  the  same  as  the  line  in  the  figure  with  a  value  of  wi  that  corresponds 
to  the  actual  value  of  tvi.  For  example,  if  the  processor  is  matched  to  a 
value  of  wi  equal  to  zero,  the  processor  output  will  be  the  wj  =  0  line  in 
the  figure.  A  very  good  estimate  of  range  is  then  obtained.  If,  however,  the 
processor  is  mismatched  to  some  other  value  of  wi,  then  a  distinctive  peak 
at  the  true  source  range  is  not  seen  in  the  corresponding  line  in  the  figure 
and  the  processor  performance  is  then  poor.  On  the  other  hand,  the  OUFP 
averages  across  the  lines  and  finds  a  line  with  a  high  peak  at  the  true  range. 
(The  line  corresponding  to  the  actual  value  of  u?i.)  This  implies  the  OUFP 
output  will  also  have  a  peak  at  the  true  source  range.  Similar  comments 
apply  to  the  depth  localization  cases. 

The  joint  range  and  W2  pdf  is  shown  in  Fig.  16.  A  ridge  appears  at  the 
true  range,  but  no  easily  identifiable  peak  appears  in  the  ridge.  This  means 
that  good  range  estimation  is  possible,  but  good  estimates  of  W2  are  not 
obtainable.  This  is  predictable  from  the  mismatched-OMFP  results.  The 
mismatched-OMFP  range  estimates  are  not  sensitive  to  uncertainty  in  loj, 
which  implies  the  joint  range  and  W2  pdf  heis  a  ridge  across  all  values  of  W2. 
This  makes  estimation  of  W2  difficult. 

The  joint  depth  and  W2  pdf  is  shown  in  Fig.  17.  There  is  a  small  ridge 
at  the  true  source  depth  but  it  is  quite  difficult  to  see  in  the  graph.  The 
low  ridge  height  indicates  difficulty  when  trying  to  estimate  depth  and  W2 
simultaneously. 

Of  course,  the  performance  is  a  function  of  the  signal-to-noise  ratio  (de¬ 
tection  factor).  Figs.  18  and  19  are  the  same  as  the  two  previous  figures 
except  the  detection  factor  has  been  raised  by  14  dB.  Since  there  is  a  peak 
in  the  joint  range  or  depth  and  W2  pdf  at  this  signal-to-noise  ratio,  good 
estimates  should  be  obtained.  Since  the  mismatched-OMFP  corresponds  to 
a  single  line  in  these  figures,  this  processor  will  be  sensitive  to  uncertainty  in 
W2  at  this  signal-to-noise  ratio. 
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5.4  Summary 

The  degradation  in  performance  that  can  occur  when  a  matched  field  pro¬ 
cessor  is  used  in  an  environment  for  which  it  has  not  been  designed  has  been 
illustrated.  In  an  uncertain  environment  the  performance  of  the  matched 
field  processor  may  be  extremely  poor.  On  the  other  hand,  the  optimum 
uncertain  field  processor  is  robust  in  an  uncertain  environment.  A  properly 
designed  uncertain  field  processor  can  achieve  good  performance  regardless 
of  the  actual  acoustic  environment. 
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6  Results-Sensitivity  to  Small  Scale  Range 
Independent  Perturbations 

In  the  previous  section  we  showed  that  the  OUFP  is  less  sensitive  to  uncer¬ 
tainty  in  the  major  parameters  of  the  SVP  than  a  matched  field  processor. 
The  assumed  profile  uncertainty  is  not  actually  typical  of  that  of  the  real 
ocean.  In  this  section  sensitivity  with  respect  to  perturbations  in  the  SVP 
on  the  order  of  those  which  might  occur  in  the  surface  layer  is  investigated 
by  Richardson  [10,  11]. 

6.1  The  Sound  Velocity  Profile  Model 

The  sound  velocity  profile  is  modeled  as 

c{z)  =  CMiz)  +  Cp  exp  (-2V2<^p)  (64) 

where  cm(z)  is  a  Munk  profile  with  parameters  Co  (axis  sound  speed)  equal 
to  1500  m/s  and  zq  (axis  depth)  equal  to  1000  m.  (This  is  the  typical  profile 
shown  in  the  introductory  section.)  The  second  term  represents  a  Gaussian 
perturbation  in  the  SVP  that  occurs  at  the  water  surface.  For  the  results 
presented  in  this  section,  the  perturbation  standard  deviation  width  is  20  m. 
The  localization  sensitivity  with  respect  to  the  height  of  the  perturbation, 
Cp,  is  investigated. 

Urick  [14]  shows  several  temperature  profiles  of  the  surface  layer  as  mea¬ 
sured  near  Bermuda  over  a  24  hour  period.  The  profiles  vary  by  about  5°C 
during  this  time  period.  (Cooling  occurs  at  night  with  warming  taking  place 
during  the  day.)  This  corresponds  to  a  variation  in  sound  speed  of  approxi¬ 
mately  8  m/s.  We  investigation  variations  of  this  order  of  magnitude. 

6.2  Estimation  Sensitivity 

As  in  the  previous  section,  the  range  to  the  source  is  50000  meters.  A  10 
element  array,  with  elements  200  m  apart,  is  used.  The  detection  factor  is 
86  dB.  The  normal  mode  acoustic  model  is  used  to  calculate  the  observation 
and  the  replica  fields. 

A  shallower  source  depth  is  used  and  we  look  at  localization  performance 
as  a  function  of  source  depth  and  frequency. 
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6.2.1  Range  Estimation 

For  the  first  set  of  results  the  source  is  at  a  depth  of  15  m  and  is  operating  at 
a  frequency  of  15  Hz.  The  mismatched-OMFP  range  estimation  sensitivity 
to  uncertainty  in  Cp  is  shown  in  Fig.  20.  {The  Cp  axis  is  labeled  Delta  in  the 
figure.)  The  mismatched-OMFP  is  matched  to  a  perturbation  height  of  zero. 
The  processor  shows  very  little  sensitivity  to  perturbation  height  uncertainty 
over  the  range  of  mismatch  examined. 

The  corresponding  OUFP  output  is  shown  in  Fig.  21.  This  OUFP  as¬ 
sumes  uncertainty  in  the  surface  layer  and  that  the  perturbation  amplitude 
is  uniformly  distributed  between  -10  and  10  m/s.  The  high  ridge  indicates 
that  this  processor  is  also  not  sensitive  to  uncertainty  in  the  perturbation 
amplitude. 

The  insensitivity  is  due  to  the  large  ratio  of  wavelength  to  perturbation 
width.  At  15  Hz,  the  acoustic  wavelength  is  approximately  100  m.  This  is  five 
times  larger  than  the  standard  deviation  width  of  the  perturbation.  Figs.  22 
and  23  show  the  mismatched-OMFP  and  OUFP  output  for  an  acoustic  source 
operating  at  60  Hz.  The  wavelength  at  this  frequency  is  approximately  25  m, 
which  is  about  equal  to  the  perturbation  width.  At  this  higher  frequency 
the  mismatched-OMFP  is  quite  sensitive  to  mismatch  greater  than  about 
2.5  m/s.  On  the  other  hand,  the  OUFP  performs  robustly  with  respect  to 
uncertainty  in  the  actual  value  of  Cp. 

The  sensitivity  of  matched  field  processing  is  a  function  not  only  of  fre¬ 
quency,  but  source  depth  as  well.  Figs.  24  and  25  correspond  to  Fig.  22 
except  the  source  depths  are  30  and  60  m.  At  30  m  the  mismatched-OMFP 
is  still  somewhat  sensitive  to  Cp  uncertainty.  At  60  m  the  mismatched-OMFP 
output  is  relatively  insensitive  to  uncertainty  in  Cp.  At  this  depth  the  source 
is  well  below  the  surface  layer  perturbation  and  the  resulting  acoustic  field 
is  relatively  unaffected  by  it. 

6.2.2  Depth  Estimation 

We  now  look  at  depth  estimation  performance  when  the  range  to  the  source 
is  known  (50000  m).  The  source  is  at  a  depth  of  15  m  for  all  of  the  results 
in  this  section. 

Fig.  26  shows  the  sensitivity  of  the  mismatched-OMFP  to  mismatch  in 
perturbation  height  for  the  15  Hertz  source.  The  peak  is  in  approximately 
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the  right  location  when  the  processor  is  matched,  but  the  resolution  is  poor. 
As  the  value  of  Cp  moves  farther  from  zero  (i.e.,  the  mismatch  between  the 
designed  environment  and  the  actual  environment  becomes  greater)  the  peak 
in  the  pdf  shifts  downward  from  the  true  source  depth.  Interestingly,  the 
resolution  increases. 

Fig.  27  shows  the  corresponding  output  of  the  OUFP.  The  peak  is  in  the 
right  location  for  all  values  of  the  perturbation  amplitude,  but  the  resolution 
is  poor. 

Figs.  28  and  29  show  the  processor  outputs  after  increasing  the  source 
frequency  to  60  Hz.  The  mismatched-OMFP  does  a  good  job  of  determining 
the  source  depth  for  values  of  Cp  near  zero  and  for  positive  values  of  Cp,  but 
the  performance  of  the  OUFP  is  clearly  superior  to  that  of  the  mismatched- 
OMFP  for  almost  all  values  of  Cp.  The  resolution  of  both  processors  is  in¬ 
creased,  due  to  the  higher  source  frequency. 

6.3  Summary 

The  performance  of  the  mismatched-OMFP  and  the  OUFP  is  dependent 
upon  source  position,  the  acoustic  environment,  and  the  frequency.  The  ro¬ 
bust  performance  of  the  OUFP,  relative  to  that  of  the  mismatched-OMFP, 
was  illustrated  in  several  cases  for  realistic  SVP  pertrubations.  At  high  fre¬ 
quencies  the  results  indicate  that  mismatched-OMFP  may  be  very  sensitive 
to  even  small  uncertainties  in  the  SVP.  The  mismatched-OMFP  is  limited  to 
applications  in  which  the  size  of  any  uncertain  SVP  perturbations  is  a  small 
fraction  of  the  acoustic  wavelength. 
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7  Results- Sensitivity  to  Range  Dependent 

Pert  urbat  ions 

In  this  section  we  investigate  and  compare  the  sensitivity  of  matched  field  and 
uncertain  field  processing  in  a  range-dependent  acoustic  environment  [10,  11]. 


7.1  The  Sound  Velocity  Profile  Model 


The  range-dependent  sound  velocity  profile  is  modeled  as 

ciz)  =  cm{2)  -I-  Cpexp  [(r  -  rp)^/(T^  +  {2-  Zp)V<^z])  (65) 


where  cm(2)  is  a  Munk  profile  with  parameters  cq  (axis  sound  speed)  equal  to 
1500  m/s  and  zo  (axis  depth)  equal  to  1000  m.  (This  is  the  same  profile  used 
in  the  previous  sections.)  The  second  term  represents  a  range-dependent 
Gaussian  perturbation  in  the  SVP  with  a  peak  at  range  Tp  and  depth  Zp. 
The  standard  deviation  width  in  range  is  cr^  while  that  in  depth  is  cTj.  The 
sensitivity  with  respect  to  Cp  is  again  investigated 

The  perturbation  peak  is  placed  at  a  range  of  25000  meters  from  the 
source  and  at  axis  depth  ( 1000  meters).  The  range  standard  deviation  width 
is  5000  meters  while  the  depth  width  is  250  meters.  The  height  of  the  peak 
is  varied  between  -5  and  5  meters/second. 


7.2  Suboptimum  Localization 

The  standard  normal-mode  model  is  not  range-dependent,  so  the  ray  theory 
model  is  used  to  calculate  the  replica  fields  and  the  observation. 

Calculation  of  the  replica  fields  for  the  optimum  uncertain  field  processor 
requires  that  a  new  ray  trace  be  performed  for  each  possible  source  position 
and  each  environment.  The  computer  processing  time  is  prohibitive  for  any 
problems  of  interest. 

Instead  of  examining  the  performance  of  the  range-dependent  optimum 
uncertain  field  processor  we  look  at  three  suboptimum  range-independent 
processors.  These  suboptimum  processors  use  range-independent  ray  trac¬ 
ing  to  calculate  the  replica  fields.  (The  principle  of  acoustic  reciprocity  can 
be  applied  in  a  range  independent  environment  to  greatly  reduce  the  num¬ 
ber  of  ray  traces  that  have  to  be  done  to  compute  the  replica  fields.  Replica 
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fields  corresponding  to  the  same  source  depth  and  environment,  but  different 
ranges,  can  be  computed  from  the  same  ray  trace.)  Since  the  actual  environ¬ 
ment  is  range-dependent,  none  of  the  suboptimum  processors  can  calculate 
a  replica  field  which  exactly  matches  the  observed  field. 

The  first  suboptimum  processor  is  the  mismatched  optimum  matched 
field  processor  (mismatched-OMFP).  This  processor  is  matched  to  the  range 
independent  environment  which  exists  in  the  absence  of  the  perturbation. 

The  second  suboptimum  processor  is  the  mismatched  optimum  uncer¬ 
tain  field  processor  (mismatched-OUFP).  It  calculates  replica  fields  from 
range- independent  profiles  which  have  a  Gaussian  perturbation  at  the  axis 
depth.  The  height  of  this  perturbation  is  cissumed  unknown  and  uniformly 
distributed  between  -2  and  2  meters/second.  Since  the  actual  range  averaged 
profile  lies  within  this  region  of  uncertainty  it  is  hoped  that  good  performance 
can  be  achieved  with  this  processor.  (If  the  actual  perturbation  height  is 
5  meters/second,  the  range  averaged  profile  (average  out  to  75000  meters) 
is  a  Munk  profile  with  Gaussian  perturbation  of  height  0.8  meters/second.) 
Note  that  this  processor  is  an  uncertain  field  processor,  but  it  doesn’t  use 
the  correct  uncertain  environment. 

The  third  suboptimum  processor  is  a  new  uncertain  field  processor.  It 
averages  conventional  ambiguity  functions  which  result  from  different  replica 
fields.  The  average  is  calculated  over  the  same  range-independent  profile 
uncertainty  as  used  by  the  mismatched-OUFP.  This  processor  is  termed  mis¬ 
matched  supoptimum  uncertain  field  processor  (mismatched-SUFP). 


7.3  Estimation  Sensitivity 

The  same  10  element  array  used  in  previous  sections  is  used  here.  The  source 
is  on  axis  (1000  m)  and  at  a  range  of  75000  m.  The  frequency  is  100  Hertz. 
The  detection  factor  is  93  dB. 

7.3.1  Range  Estimation 

We  consider  the  problem  of  estimating  only  the  range  coordinate  when  the 
source  depth  is  known. 

The  mismatched-OMFP  output  as  a  function  of  actual  perturbation  height 
is  shown  in  Fig.  30.  Note  that  when  Cp  (labelled  Delta  in  the  figure)  is  zero 
the  actual  environment  is  range-independent  and  the  mismatched-OMFP  is 
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Figure  30:  Mismatched-OMFP  a  posteriori  range  pdf  output  as  a  function 
of  perturbation  amplitude.  The  true  source  range  is  75000  m. 


matched  to  the  environment.  In  this  region  strong  peaks  near  the  true  source 
range  of  75000  m  can  be  observed.  When  the  absolute  value  of  the  actual 
perturbation  height  is  greater  than  about  2  m/s  ambiguities  appear  in  the 
output  and  localization  performance  decreases  rapidly. 

The  output  of  the  mismatched-OUFP  is  shown  in  Fig.  31.  This  processor 
performs  slightly  better  than  the  previous  one.  Some  of  the  large  spurious 
peaks  have  been  suppressed  and  the  ridge  at  75000  m  is  slightly  wider. 

The  output  of  the  mismatched-SUFP  is  shown  in  Fig.  32.  The  perfor¬ 
mance  of  this  processor  is  quite  poor  for  almost  all  values  of  the  pertur¬ 
bation  amplitude.  The  figure  illustrates  that  seemingly  reasonable  ad  hoc 
approaches  to  the  design  of  localization  algorithms  may  result  in  processors 
with  very  poor  performance. 

7.3.2  Depth  Estimation 

We  now  consider  the  problem  of  estimating  the  depth  coordinate  when  the 
range  is  known. 

The  output  of  the  mismatched-OMFP  as  function  of  perturbation  height 
is  shown  in  Fig.  33.  There  is  a  peak  at  the  true  source  depth  when  the 
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Figure  33:  Mismatched-OMFP  a  posteriori  depth  pdf  output  as  a  function 
of  perturbation  amplitude.  The  true  source  depth  is  1000  m. 

perturbation  amplitude  is  zero  (no  mismatch),  but  the  peak  of  the  pdf  shifts 
away  from  the  true  depth  as  the  perturbation  amplitude  moves  away  from 
zero.  The  peak  in  the  pdf  appears  at  a  depth  almost  200  m  away  from  the 
true  depth  ath  the  extreme  values  of  the  perturbation  amplitude. 

The  corresponding  output  of  the  mismatched-OUFP  is  shown  in  Fig.  34. 
The  performance  improvement  over  the  mismatched-OMFP  case  is  more  no¬ 
ticeable  here  than  in  the  range  estimation  case.  Some  spurious  peaks  are 
present  in  the  figure,  but  there  is  a  ridge  very  near  the  true  source  depth 
across  all  values  of  the  perturbation  amplitude. 

The  mismatched-SUFP  output  is  shown  in  Fig.  35.  This  processor  per¬ 
forms  much  better,  relative  to  the  mismatched-OMFP,  than  in  the  range 
estimation  Ccise.  A  ridge  appears  at  the  correct  depth  over  the  portion  of 
the  figure  near  where  the  perturbation  amplitude  is  zero,  but  large  spurious 
peaks  appear  away  from  the  true  source  depth  as  the  perturbation  amplitude 
increases. 

To  put  the  performance  of  these  suboptimum  processors  in  perspective, 
the  output  of  the  optimum  uncertain  field  processor  for  this  environment  is 
shown  in  Fig.  36.  This  processor  assumes  the  environment  is  range-dependent 
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Figure  36;  OUFP  a  posteriori  depth  pdf  output  as  a  function  of  perturbation 
amplitude.  The  true  source  depth  is  1000  m. 


with  a  Gaussian  perturbation  of  unknown  amplitude.  The  range- dependent 
acoustic  ray  propagation  model  is  used  to  calculate  the  replica  fields.  The 
optimum  processor  is  quite  robust  with  respect  to  uncertainty  in  the  pertur¬ 
bation  height. 

7.4  Summary 

In  this  section  we’ve  compared  the  performance  of  three  suboptimum  pro¬ 
cessors.  The  three  processors  used  a  range-independent  propagation  model 
to  calculate  the  replica  fields,  although  the  actual  environment  is  range- 
dependent.  The  suboptimum  processor  which  is  beised  on  the  optimum  un¬ 
certain  field  processor  had  the  best  performance.  Although  this  processor 
doesn’t  perform  as  well  as  the  optimum,  the  reduced  processing  time  require¬ 
ments  make  it  an  attractive  alternative  to  the  optimum  in  some  situations. 
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